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We derive exact formulas for bipartite von Neumann entanglement entropy after partial projective 
local measurement in 1 + 1 dimensional conformal field theories with periodic and open boundary 
conditions. After defining the set up we will check numerically the validity of our results in the case 
of Klein-Gordon held theory (coupled harmonic oscillators) and spin-1/2 XX chain in a magnetic 
held. The agreement between analytical results and the numerical calculations is very good. We 
also hnd a lower bound for localizable entanglement in coupled harmonic oscillators. 
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In the last couple of decades bipartite entanglement 
entropy attracted a lot of attention in the high energy 
physics and the many body condensed matter physics 
mostly due to the area law property of bipartite von Neu¬ 
mann entanglement entropy which is the most celebrated 
measure of entanglement [l|, l2| . In one dimensional quan¬ 
tum systems while the area law is usually valid just for 
massive (gapped) systems for the critical systems the en¬ 
tanglement entropy of subsystem with the size I of the 
ground state follows the logarithmic formula S = ^ log/, 
where c is the central charge of the underlying conformal 
field theory (CFT) describing the low energy behavior of 
the critical system Since the central charge of the 

system usually can fix the universality class of the sys¬ 
tem calculating entanglement entropy for a system can 
give a lot of insight about the possible universality class 
of the system. In particular, since numerical calculation 
of the entanglement entropy by using the techniques of 
DMRG in one dimension is now a well-known method @ 
one can easily study the critical and non-critical proper¬ 
ties of the quantum systems in one dimension by study¬ 
ing the entanglement entropy. It is worth mentioning 
that entanglement entropy is not the only quantity which 
gives directly the central charge of the system. AmoM 
other measures one can name, mutual informations [7H9| 
and quantum contours 0 , see also other related works 
[III Il2j | . The former one is based on local measurements 
in particular basis, so called conformal basis 0 , 0 . One of 
the simple properties of the entanglement entropy which 
makes it more appealing in numerical calculations with 
respect to other measures such as Shannon and Renyi 
mutual informations is that it is completely independent 
of the basis that we write the wave function. However, 
since it is a very non-local quantity it is not an easy quan¬ 
tity to be measured in the experiment, for recent devel¬ 
opments see (l3l - [i0 . Although bipartite entanglement 
entropy has been studied in length in many quantum 
systems there are few studies regarding entanglement in 
multipartite systems m. One of a few entanglement 
measures regarding tripartite systems is localizable en¬ 
tanglement introduced in [I^, see 0 for review. The lo¬ 
calizable entanglement is defined as the maximal amount 
of entanglement that can be localized, on average, by do¬ 


ing local measurements in part of the system. In other 
words the localizable entanglement between the two parts 
B and B after doing local measurement in the rest of the 
system A is defined as 

EiociB,B) = sups^PkE{\'4)k)bb) ( 1 ) 

K 

where S is the set of all possible outcomes 
((pk, £’(|^if)g^) of the measurements and E is 
the chosen entanglement measure. Because of the 
maximization over all the possible local measurements 
the localizable entanglement is a very difficult quantity 
to calculate. However, in those cases that B and B are 
single spins and E{.) is the concurrence one can derive 
interesting lower bounds to the localizable entanglement, 
see 0 and references therein and for other related works 
see also M- The quantity has been also measured 
experimentally in a system of two photons in a noisy 
surrounding in [^ . Due to the complexity of the 
definition of localizable entanglement, to the best of 
our knowledge, the case of B and B not being just 
single particles but many particle subsystems has not 
been addressed in the literature. In this work instead of 
calculating the very complicated quantity of localizable 
entanglement we calculate E{\'ipK) bb) appearing in the 
definition ([S]) for particular natural basis. In the case of 
harmonic oscillators our exact results are useful to find 
lower bound for the localizable entanglement. 

Consider the ground state of a many body system (for 
example a spin chain) written in particular basis, i.e. 
\ipg) = where for example in the case of spin 

chain I are all the possible spin configurations in the 
direction. If we do local projective measurement of a local 
quantity (for example in the spin chain) in some part 
of the system that part will take definite spin direction 
and will be disentangled from the rest. In other words 
one can write IV'g) = J^i^iMkII) = J2jO'Jk\J)\K), 
where Mk is the projective measurement operator of the 
subsystem with the outcome K and the sum is now over 
the spins of the rest of the system. Notice that if we start 
with the ground state of a spin chain the resulting state 
after measurement is nothing to do with the ground state 
of the rest of the system. The entanglement content of 
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the new wave function will be dependent on the basis 
that we chose for our measurement and of course to the 
outcome of the measurement, in other words, the entan¬ 
glement content of the |'0g) is dependent on Mk- In a re¬ 
cent development 0 it is shown that although, in general, 
doing measurement in arbitrary basis is not compatible 
with CFT set up there are some natural basis (confor¬ 
mal basis) that are closely related to boundary CFT. In 
other words if one makes a measurement in these basis 
one can still hope to preserve some of the universal prop¬ 
erties of the system and be able to calculate the entan¬ 
glement entropy. We noticed that engineering quantum 
spin chains and making projective measurement on nat¬ 
ural local basis is now possible with optical lattices and 
ion trapping techniques, see [ 2 l| and references therein. 
However, calculating entanglement entropy after partial 
local measurement in many body quantum chains has 
not been investigated yet. Analytical calculation of bi¬ 
partite von Neumann entanglement entropy after projec¬ 
tive measurement in natural conformal basis is the main 
purpose of this article. To do that we first define our set 
up and fix our assumptions and then we will analytically 
calculate entanglement entropy using twist operators of 
CFT. Exact formulas will be derived and then we will 
check the validity of our results by explicit examples in 
the field theories such as Klein-Gordon field theory (cou¬ 
pled harmonic oscillators) and in quantum spin chains 
such as, XX model in a magnetic field. 


I. BIPARTITE ENTANGLEMENT ENTROPY 
AFTER PARTIAL MEASUREMENT IN CFT. 


The basic setup for our problem is as follows: con¬ 
sider the ground state of a periodic critical system in one 
dimension with the total size L. Then we make a pro¬ 
jective local measurement in the subsystem A with the 
length s of this system in particular local basis. After the 
measurement the subsystem A will be decoupled (unen¬ 
tangled) from the complement A. However, the wave 
function of the subsystem A after measurement, i.e. 
is highly entangled. In other words if we take a subsys¬ 
tem B C Ait will be entangled with B, where BUB = A. 
From now on we will take the length of the subsystem B 
equal to I and for simplicity we will take B and H in a way 
that they are adjacent simply connected domains. Here 
we are interested to calculate the entanglement entropy 
of the subsystem B with respect to B, i.e. Sb, in a criti¬ 
cal quantum chain. To calculate such a quantity it is first 
useful to think about the Euclidean version of the system 
after measurement. For a periodic boundary conditions 
using the transfer matrix approach one can simply map 
the system after measurement to a cylinder with a slit on 
it as it is demonstrated in Figure 1, see Q- In general 
this system is not necessarily an example of boundary 
CFT except for those cases that the subsystem A after 
measurement picks up a configuration which renormal¬ 
izes to a boundary CFT. We will discuss this important 



FIG. 1: Euclidean version of the quantum chain with total 
length L. The removed slit domain A has length s and we are 
interested in calculating entanglement entropy of the region 
B with length s with the complement in the quantum spin 
chain. The twist operator can be put at point P. 


issue later in the case of critical spin chains. Now if we 
consider that the outcome of the measurement leads us 
to a boundary CFT, then based on the Cardy-Calabrese 
argument @ the calculation of the entanglement entropy 
Sb should boil down to the calculation of the one point 
correlation function of twist operator sitting on the bor¬ 
der between B and B on the cylinder. In other words we 
have: 

Pb=- lim ^ < %{P) (2) 


where TP. is the twist operator with conformal weight 
hn = if (1 ~ fr) and c is the central charge. The one 
point function of the twist operator on the cylinder with 
a slit can be easily found by mapping the system to the 


upper half plane by using the map z{w) = 

and the following well-known formula of CFT f^ : 


< >cyl/sHt~ i^wZ) < Tn >H, (3) 


/ 

where < Tn >h= ( 2z(s/2-i-;)) ) point func¬ 

tion on the upper half plane %. a is the lattice spacing 
and we took the coordinate of the twist operator at § -I-/. 
Putting all together we will have 


SB = ^ln(- 
6 V TT 


L sin y{ 1 + s) sin fl ' 


+ 7i 


(4) 


where 71 is a constant and the dots are subleading 
terms. Notice that in the above formula s can not go to 
zero because cylinder with a slit is topologically different 
from cylinder. To get the result for before measurement 
case one can simply put s = a which is the smallest scale 
in the system, then we will have Sb = f In-(^sin^) 
which is the well-known result for the bipartite entan¬ 
glement entropy in CFT Q . In the next sections we will 
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verify the validity of the above equation in two important 
cases: Klein-Gordon field theory (harmonic oscillators) 
and the critical XX model. 


II. HARMONIC OSCILLATORS. 

In this section we would like to calculate bipartite en¬ 
tanglement entropy of the ground state of coupled har¬ 
monic oscillators after measuring the position of a string 
of oscillators. Consider the Hamiltonian of L-coupled 
harmonic oscillators, with coordinates (pi,... ,(j)L and 
conjugated momenta tti, ..., ttl: 



n—1 


L 

E ' 

n,n' — l 


(5) 


The ground state of the above Hamiltonian has the fol¬ 
lowing form 


^0 = ( 




( 6 ) 


FIG. 2: Entanglement entropy of subregion B for short-range 
coupled harmonic oscillators with total length L = 50 and the 
measurement region sizes s = 10,16 and 20 with respect to 
ln/(L,s,Z), where f{L,s,l) = § ■ The full line 

is the function with c = 1 and 71 = 0 . 21 . 


to find a lower bound for the localizable entanglement in 
this system as follows 


One can calculate the two point correlators Xjj = 
tr {p£)(pipj) and Pp = tr using the K matrix, 

where pu is the reduced density matrix of domain D. 
The squared root of this matrix, as well as its inverse, 
can be split up into coordinates of the subsystems D and 
D , i. e.. 




( Xd Xob\ 


^1/2 



Pdd\ 

Pd)- 


Sioc{B,B) > - In 
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1 /L sin + s) sin ^ 


+ 71- ( 8 ) 


Notice that the above result is based on this fact that 
the scaling properties of the entanglement entropy after 
partial measurement in harmonic oscillators is indepen¬ 
dent of the outcome of the measurement. This is not 
necessarily correct in generic systems. 


The spectra of the matrix 2C = \/Xr>Pn, can be used to 
calculate the entanglement entropy, see and reference 
therein, 


S' = tr 


(C+i)log(C + i)-(C'-i)log(C-i) 


(7) 


Now if we do measurement on the position of all the os¬ 
cillators {(pi] G A they will take some definite values (for 
example, {pi}=0 ) and eventually will get decoupled from 
the rest of the oscillators. In other words the final state 
will be the same as m but instead of we need to 
consider which is a subblock of the matrix 

corresponding to the oscillators in the subsystem A. This 
means that we now have a new Gaussian wave function 
and one can calculate its bipartite entanglement entropy 
with the formula (l7|). The results for the short-range har¬ 
monic oscillators (discrete Klein-Gordon field theory) are 
shown in the figure 2 which show perfect agreement with 
the equation (|4]) if we consider c = 1 which is the central 
charge of the free field theory. Since the above results are 
independent of the outcome of the measurement of the 
the position of the oscillators {(pi) € A one can use them 



FIG. 3: Entanglement entropy of subregion B with length I 
after measurement on subsystem with length s with ferromag¬ 
netic outcome for XX model in a magnetic field . From top to 
bottom we have n/ = |, i and i. In our numerics Z -|- s = 40 is 
fixed. The full line is the function with L —^ 00 and c = 1 
and 71 = 0.53 and the arrows are at = (1 — n/)(Z - 1 - s) 
from left to right for n/ = |^ and | . 
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III. XX MODEL IN A MAGNETIC FIELD. 

In this section we would like to calculate bipartite en¬ 
tanglement entropy after partial projective measurement 
in XX model in a magnetic field. The method that we 
use can work equally for all the other spin chains that 
can be mapped to free fermions. The Hamiltonian of XX 
model is as follows: 

H = -\Y^ (o-fcrf+i -b erf crf+i - 2/lcrf) . (9) 

After using simple Jordan-Wigner transformation, i.e. 
Cl = i Hamiltonian will have the fol¬ 

lowing form 

H = -Y^ (cjci+i + + 2h{cjci - i)) . (10) 

The entanglement entropy of a subsystem can be calcu¬ 
lated easily, see , by first calculating the reduced den¬ 
sity matrix pB = Ae” ^and then diagonalizing 
it, the final formula is 



FIG. 4: Entanglement entropy of subregion B with length 
I after measurement on subsystem with length s with anti¬ 
ferromagnetic outcome for XX model in half filling . In our 
numerics Z -b s = 40 is fixed. The full line is the function @ 
with L —>■ oo and c = 1 and 71 = 0.61. 
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-tr 


(l-C) log(l-C')-bClogC 


( 11 ) 


The matrix H is related to the correlation matrix C 
a.s H = log(C'“^ ~ !)■ In our case the correlation ma¬ 
trix for infinite system has the form Cij =< c|cj >= 
! where nj = Looking to the 

Jordan-Wigner transformation it is easy to see that any 
measurement in the basis on particular site can be 
translated to the measuring the number of fermions in 
that site. In other words if one measures the cr^ in all 
the sites of a subsystem with size s the outcome will 
be one of 2^ possible configurations which can be eas¬ 
ily translated to the conhgurations made of presence or 


absence of fermions in the sites of the subsystem. For 
simplicity we first consider that the outcome of the mea¬ 
surement is a string of s fermion occupied sites. This 
can be calculated easily by using Grassmann numbers 
by first calculating the reduced density matrix for a sub¬ 
system with length I + s and then finding the reduced 
density matrix of the subsystem with length I with the 
assumption that the outcome of the measurement in the 
subsystem s is a string of filled sites. The method is ex¬ 
plained with full detail in the appendix. The results for 
different filling factors are shown in the Figure 3. There 
are some comments in order: first of all because of the 
C/(l) symmetry of the XX model the number of particles 
in the system is conserved and for this reason as far as 
n/ is small it is very difficult to have a string of sites with 
fermions. In other words for example in the XX model 
with h = 0 the most probable outcome is an antiferro¬ 
magnetic string rather than ferromagnetic string. In Q 
it is already conjectured that for the half filling case most 
probably the ferromagnetic configuration is not going to 
lead to a boundary conformal field theory and so the very 
first assumption that we used is going to fail. Despite this 
argument we found surprisingly that for I > (A— I)s the 

formula (|4]) works perfectly. 

One can also check the results for those cases that the 
outcome of the measurement is antiferromagnetic string. 
It is expected that this case leads to Dirichlet boundary 
condition in the bosonization language and so it is related 
the boundary CFT. The results presented in the Figure 
4 indeed show that the formula @ works perfectly also 
in this case. Since the ferromagnetic and antiferromag¬ 
netic configurations are at the two extreme sides of the 
all possible configurations based on the above numerical 
results one can conjecture that independent of the out¬ 
come of the measurement in the basis the formula (l4|) 
will work if the I is bigger enough than s. Note that Since 
the number of conhgurations increases exponentially this 
conjecture is beyond what we can check numerically. 


IV. OPEN BOUNDARY CONDITIONS 


One can also do all the CFT calculations in the case 
of open boundary condition, the only difference is that 
now we have a strip with a slit instead of a cylinder. One 
can simply map a strip with a slit to upper half plane by 

using the map z{w) = y^I — ^2 M ■ After some algebra 
we have 


^B=gln( 


/ 2L cos ^ - cos TT^ 7t{1 + s) 


a cos 


2 

2L 


cot • 


2L 


72.(12) 


The limit s —>■ 0 simply gives us the well-known formula 
Sb = I Insin ^), see We have checked the validity 
of this formula in the case of short-range harmonic oscilla¬ 
tors using the same method that we hired for the periodic 
BC. The results shown in Figure 5 are fairly compatible 
with the formula m- We noticed that although the log 
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formula is perfectly compatible there is almost 8 percent 
deviation from central charge that might be due to finite 
size effects. 



FIG. 5: Entanglement entropy of subregion B for a sys¬ 
tem (short-range coupled harmonic oscillators) with total 
length L = 50 and the measurement region sizes s = 
0,10 and 20 with respect to In f{L, s,l), where f{L,s,l) — 

— cot . The full line is the function (IT^ 

TT a r.os-^ 4^ 2L ' 


with c = 1 and 72 = — 0 . 02 . 


V. CONCLUSIONS. 

We described a general set up for calculating bipartite 
entanglement entropy after local projective measurement 
in critical one dimensional quantum systems. Exact for¬ 
mulas were derived for bipartite entanglement entropy af¬ 
ter “conformal measurements “ in the case of periodic and 
open boundary conditions. The formulas were checked in 
explicit examples of free bosonic system and XX model 
in a magnetic field. We noticed that since bosonization 
of XXZ model in the basis leads to free bosonic sys¬ 
tem the antiferromagnetic outcome of the measurement 
in this basis should be compatible with the results of free 
bosonic system Q. In the case of XX model in a mag¬ 
netic field we showed that if we do our measurements 
in the cr^ basis independent of having ferromagnetic or 
antiferromagnetic outcome for our measurements the bi¬ 
partite entanglement entropy in particular regimes can 
be described with GET formulas. We also derived a 
lower bound for the localizable entanglement in the case 
of harmonic oscillators. There are interesting questions 
remained to be answered: first of all in our free fermion 
approach we were able to handle just basis, it is im¬ 
portant to look to the other basis especially cr^ basis by 
using exact diagonalization methods. Another interest¬ 
ing question is investigating the same questions in the 
case of non-critical quantum systems and the validity of 
the area law. Some of these questions will be discussed 
in a forthcoming paper [^ . Finally understanding the 
problem in the holographic set up (26| will surely help to 
extend some of these results to higher dimensions. 
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Appendix 


Bipartite entanglement entropy after partial 
measurement of fermion occupation numbers in free 
fermions 


In this supplementary note we would like to present a 
method which can be used to measure the entanglement 
entropy after measuring the fermion occupation numbers 
in a subsystem in a generic free fermion system. The 
method is based on Grassmann numbers and it is a gen¬ 
eralization of the work Although the method can 

be also used for the most generic free fermion system [ 2 ^ 
here we will concentrate on a system with the following 
reduced density matrix: 

PS = (SI) 

The matrix H is related to the correlation matrix C as 
H — log(C'“^ — 1). Finally one can find the entanglement 
entropy using the following formula, see for example ; 


Sb 


-tr 


(l-C) log(l-C')-f ClogC 


(S2) 


To calculate the reduced density matrix after measure¬ 
ment we need to first define fermionic coherent states. It 
can be defined as follows 


1C >= ICiC2...Civ >= >, (S3) 

where Ci’s are Grassmann numbers following the proper¬ 
ties: CnCm + CmCn = 0 and = 0. Then it is easy 

to show that 


g|C>=-C.IC>. (S4) 

Using the above formula one can derive the following use¬ 
ful formula: 

< (S5) 

With the same method one can also define another kind 
of fermionic coherent state as 


\rj >= \Tjir] 2 -VN >= e (S6) 

where ry^’s are Grassmann numbers. Then it is easy to 
show that 


4\v>=-Vi\n> (S7) 


and consequently 

< .^|gEy 


(S8) 
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The two formulas (IS5I) and (IS 8 I) are the main formu¬ 
las one can use to calculate reduced density matrix after 
measurement. We first discuss how the method works 
for a string of empty sites. A string of empty sites means 
that in the equation (IS5I) we need to put all the ^i’s with 
i inside the measuring subsystem equal to zero. This is 
equivalent to keeping just all the elements of the matrix 
outside the subsystem. Then we take the logarithm 
of this matrix and find a new matrix A which is going 


to be our new entanglement Hamiltonian and then one 
can easily calculate the entanglement entropy of that sys¬ 
tem with the method we explained at the beginning of 
this section. To do the same calculation for an string 
of occupied sites instead of equation (ISSl) one needs to 
use equation (ISSl) and follow the same method as we just 
explained. In principle this method can be applied for 
any possible configuration by just appropriate using the 
equations (ISSl) and (IS 8 |) . 
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